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ABSTRACT 



Fluctuations of the gamma-ray background (GRB) 
are investigated. It is assumed that the GRB is pro- 
duced by discrete sources. Two basic classes of ob- 
jects as dominant contributors to the GRB are con- 
sidered: a) blazars and b) starburst galaxies. Pre- 
dicted log N — log S counts are calculated for these 
two models. In both cases strong cosmic evolution is 
needed to produce a substantial fraction of the ob- 
served GRB. Since on the average blazars are several 
orders of magnitude more luminous than starburst 
galaxies, their number required to saturate the GRB 
is much smaller than the number of starburst galax- 
ies. In effect, fluctuations of the GRB in the blazar 
model are relatively high. It is shown that moder- 
ate improvement in the measurement accuracy of the 
GRB flux will allow us to discriminate between these 
models. 

Key words: background radiation, gamma-ray 
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1. INTRODUCTION 



Using the present detection techniques we are un- 
able to define a true nature of the diffuse extragalac- 
tic gamma-ray background. The status of the GRB 
is now analogous to that of the X-ray backgrund in 
1960s and 70s. Although some classes of objects are 
known 7-ray emitters, it is by no means certain what 
fraction of the GRB is produced by any particular 
classes of objects and even by the discrete sources in 
general. Furthermore, it seems likely that some frac- 
tion of the diffuse emission considered previously to 
be of the extra galactic origin is pr o duced in the halo 
of our Galaxy ( Btrong ct al., 1999 ). Par fe Dc Rujuk 
postulate that most if not all of the GRB is 
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produced in the Galactic halo. 

One should expect that structural features of the 
GRB are related to its origin. It is likely that trully 
diffuse processes produce the GRB of smooth distri- 
bution without strong well defined fluctuations. On 



the contrary, small scale variations would indicate 
the discrete nature of the background. In the present 
paper we concentrate on a class of models which as- 
sume that the whole GRB is produced by a well de- 
fined population of discrete sources. We discuss two 
cases in detail. The source counts are estimated and 
prospects for the detection of individual sources by 
the present-day and the near future instruments are 
assessed. Next we estimate amplitude of the GRB 
fluctuations in various angular scales characteristic 
to both these models. 



2. DISCRETE SOURCE MODEL 



It is assumed that the extragalactic 7-ray back- 
ground is produced by a population of unresolved 
discrete sources. We assume further that these 
sources are distributed randomly i.e. effects of 
source clustering on the GRB fluctuations are ne- 
glected (it is straightforward to show that for any 
realistic model of the source clustering the latter as- 
sumption is reasonable). Since it is not clear which 
class of sources contributes most to the GRB, we 
consider two types of objects which potentially could 
generate the GRB: a) blazars and b) starburst galax- 
ies. Blazars are know n to be 7-ray emitters (e.g. 
Mukherjee et al. 1997|) , while starburst have been 



considered as potential sources (e.g. Boltan fc Juch 
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icwicz 1999). The low sensitivity and poor angu 
lar resolution of the present-day 7-ray telescopes al- 
lows only for the detection of the most luminous ob- 
jects, mostly blazars. However, the contribution of 
blazars to the GRB remains highly uncertain and 
depends strongly on the cosmic evolution. The ques- 
tion whether the GRB could be produced by the 
starburst galaxies also cannot be answered directly. 
Locally these objects are much weaker than blazars, 
but their space density is higher and there are in- 
dications that these obj ects also are subject to the 
strong cosmic evolution (Saunders et al. 1990, Pear 



son & Rowan- Robinson 1996) 
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3. SOURCE COUNTS 



The source counts are related to the luminosity func- 
tion (LF) by the formula: 



N(>S) = 




(f)(L,z)dL, (1) 



where S and L denote the 7-ray flux and luminos- 
ity, respectively, (f>(L, z) is the evolving with redshift 
luminosity function and all the symbols have their 
ususal meaning. Low luminosty limit in the sec- 
ond intergral is defined as L min ( z S ) — AitSD 2 L (\ + 

where Dl is the luminosity distance and a 
is the average energy (E) spectral index of the 7-ray 
sources: L(E) ~ E a . In the subsequent calculations 
data taken from the literature have been scaled to the 
Hubble constant H a = 100kms _1 Mpc _1 and q = 0. 



f(>100MeV) [ph/(s cm ! )] 

-15 -14 -13 -12 -11 -10 -9 -8 -7 -6 




S(lOMeV) [MeV/(s cm 2 MeV)] 



3.1. AGN model 

Predicted log N — log S counts of the 7-ray blazars 
have been calculated using data on the local lumi- 
nosity function and t he luminosity evolution rat e 
of these objects from Chiang fc Mukherjee (1998| ). 



AGNs loud in the 7-rays detected by EGRET have 
been used. The local LF shows pronounced flatten- 
ing at luminosities below ~ 10 46 ergs _1 . Estimates 
of the luminosity evolution parametrized by a power- 
law (1 + z) 13 give /3 = 2.7 at redshifts smaller than 
2.5. At higher redshifts the evolution slows down 
and in the present calculations we assumed (3 = 0. 
Although both the LF and the evolution are sub- 
ject to large uncertainties, the overall shape of the 
log N — log:*? predicted hy the model only weakly dc- 



pen lds on the particlular choice of parameters. The 
counts shown in Fig. [j] are labelled 'AGN'. According 
to the present model counts flatten substantially at 
fluxes of ~ 10~ 8 ph/(s cm 2 ). It implies that most of 
the GRB is produced by relatively small number of 
bright sources. 



3.2. Star burst galaxy model 

Significantly different predictions concerning the 
source counts and the structure of the GRB are ob- 



tai ned in the model of the 7-ray emission by star 



burst galaxies proposed by Soltan & Juchniewicz 
K1999D . In this case the 7-rays are produced as a 
res ult of the inverse Compton scattering of the far 



infrared (FIR) photons on the cosmic ray electrons. 
Both cosmic microwave background photons and in- 
trinsic galaxy radiation contribute to the produc- 
tion of 7-rays. Since cosmic ray electrons manifest 
their presence in the galaxy via synchrotron radia- 
tion, measurements of radio luminosities of starburst 
galaxies provide data on the number of these elec- 
trons. Consequently, the correlation between the 
FIR and radio luminosities observed for the star- 



Figure 1. Counts of the j-ray sources producing 
the GRB. Energy flux of the lower abscissa axis and 
the photon flux of the upper axis are scaled assum- 
ing the average photon index of -2.1. The curve la- 
belled 'AGN' shows the logA^ — log S counts if the 
GRB is dominated by blazars; curves 'DE' and ! LE" 
represent counts predicted in the starburst galaxies 
model with the density evolution and luminosity evo- 
lution, respectively. Verical line marks typical detec- 
tion threshold of EGRET; source counts which have 
common envelope line with a slope of — 1 produce 
similar fraction of the GRB. 



for relatively safe estimates of the 7-ray luminosity 
of each galaxy and the local volume emi ssivity of 



the whole starburs t galaxy population (see Soltan 



Juchniewicz (199E) for details). In fact, the obtained 



figures are just lower limits for the high energy emis- 
sion of these objects. This is because apart from the 
inverse Compton mechanism, one should also take 
into consideration Brcmsstrahlung emitted by cos- 
mic ray electrons in the interaction with the galactic 
gas. 

The FIR evolution rate of starburst galaxies is not 
well determined, but some direct estimates at low 
redshifts are available (e.g. Saunders et al. 1990| ). 
Various observations provide data on the star forma- 
tion history covering redshifts up tu ~ 4 (e.g. Moor- 
wood 199q) . There is an agreement that star forma- 
tion rate (SFR) below z = 1 — 1.5 varies roughly as 
(l+z) 3 - 5 , and for higher z is either sta ble or smoothly 



dccaycs to low levels at z > 4 (e.g. Franceschini ct 
al. 1997, Hopkins et al. 2000j ) . Approximate data on 



burst galaxies (e.g. Chi fc Wolfcndalc 1990) allows 



SFR do not constrain strongly our estimates of the 
evolution of the 7-ray volume emissivity and source 
counts. This is because the 7-ray source counts de- 
pend strongly on the evolution type of the FIR prop- 
erties of starburst galaxies. If these objects are sub- 
ject to the density evolution (DE), the normalization 
of the FIR LF increases with redshift, but luminosi- 
ties on the average do not evolve. In the luminosity 
evolution (LE) case, the space concentration of ob- 
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jects is constant while their FIR luminosities (Lpm) 
increase with redshift. 

High energy luminosity of the individual galaxy due 
to the inverse Compton effect is proportional to the 
total number of cosmic ray ele ctrons and the energy 



density of FIR photons (e.g. Rybicki & Lightman 
1970). Since production of the cosmic ray electrons 
is also related to the L fir (indirectly via rate of su- 
pernova explosions), the 7-ray luminosity, L 7 , should 
be proportional roughly to the square of the Lfir. A 
more scrupulous analysis of the relevant correlations 
between radio luminosity, cosmic ray electron spec- 
trum and strength of the magnetic field, B, as well 
as radio and FIR luminos ities, and the correlation 
between the B and L FIR ( |Chi fe Wolfendale 1990| ) 



^ives the number of electro ns proportional to L u F f B 
( |Soltan fc Juchniewicz 199S| ) , and finally L 1 ~ L F f^- 
Thus, in the case of the LE of FIR luminosity func- 
tion, we expect a much stronger growth of the 7-ray 
emissivity by the starburst galaxies than in the DE 
case. Accordingly, the predicted counts for the LE 
model are much steeper. In Fig. |l| curves labelled 
'LE' and 'DE' show the count for Doth models (in 
both cases evolution rate of the FIR volume emissiv- 
ity is the same). 



4. FLUCTUATIONS 
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Sreekumar et al. (1998) using the EGRET data 
above 100 MeV have estimated the isotropy of the 
GRB by measuring the backgound intensity over 
the sky (far from the galactic center and the galac- 
tic plane). Analysis of 28 regions of approxi- 
mately 0.26 sr each gave the average flux of 1.36 x 
10 -5 ph (scm 2 sr) -1 with the rms scatter between 
the regions of 0.23 x 10~ 5 ph (scm 2 sr) -1 . The av- 
erage uncertainty of the measurement for each field 
was equal to 0.21 x 10 _5 ph(scm 2 sr) -1 . This ap- 
parent concordance implies that the detected fluctu- 
ations are consistent with the noise of the EGRET 
measurements. Below we estimate amplitude of the 
GRB fluctuations generated in the present models 
and compare our predictions with the parameters of 
the INTEGRAL instruments. 



distinguished directly. Since the predicted sensitivi- 
ties of the main INTEGRAL instruments, IBIS and 
SPI, are comparable to EGRET (assuming the aver- 
age source spectral slope of —2.1), it is unlikely that 
dicrete source observations would allow an unam- 
biguous differcntaition of both models. However, one 
can assess the source counts below the threshold of 
5(10 MeV) = 10~ 6 MeV/(scm 2 MeV) (what roughly 
corresponds to /(> 100 MeV) = 10~ 7 ph/(s cm 2 )) us- 
ing the fluctuation analysis. The number of sources 
per steradian producing most of the GRB in the 
AGN model is relatively small. In this case fluctu- 
ation of the GRB resulting from its discrete nature 
are substantially larger than in the starburst model. 

To illustrate a potential of the GRB fluctuation anal- 
ysis we consider the AGN model shown in Fig. |l|. 
Integrated counts in this case give the GRB flux of 
- 1.6 x 10- 4 MeV/(scm 2 MeVsr) at 10 MeV, what 
amounts approximately to 73 % of the actually ob- 
served total background. Roughly half of the model 
flux is produced by sources below the EGRET detec- 
tion threshold (assuming the energy index of —2.1). 
The relative rms fluctuations produced by the Pois- 
sonian distribution of sources in the 0.26 sr bins are 
of the order of 5 % of the total (observed) GRB. How- 
ever, if we use regions comparable with the SPI field 
of view (~ 0.06 sr) the rms amplitude amounts to 
10 % and for the IBIS fov of 0.025 sr it reaches 15 %. 
This figure is roughly equal to the signal-to-noise 
ratio reached in the EGRET measurements in the 
0.26 sr fields. Thus, if the INTEGRAL instruments 
give similar S/N, some estimates of the GRB fluctu- 
ations would become feasible. 

Numerical results in the above calculations depend 
strongly on the assumed model parameters. In par- 
ticular, if the blazars produce 100% of the GRB (not 
just 73%, as assumed in the example), the ampli- 
tude of fluctuations could be substantially larger. On 
the other hand, the starburst galaxy model provides 
more definite results on the amplitude of the GRB 
fluctuations. The variations of the GRB predicted 
for this latter model are very low and remain invari- 
ably below the detection threshold for a wide range 
of model parameters. 



Although integrated fluxes of the GRB predicted by 
the 'AGN' and 'starburst LE' models are similar, 
the source counts in these two cases are distinctly 
different. The vertical dashed line shows the thresh- 
old sensitivity of EGRET (~ 10~ 7 ph/(s cm 2 ) above 
100 MeV). A number of sources above this thresh- 
old have been in fact detcted by EGRET. However, 
these sources produce a small fraction of the total 
GRB. In the 'starburst LE' model most of the GRB 
is produced by sources fainter by several orders of 
magnitude, while in the 'AGN' model the substan- 
tial fraction of the GRB is generated by sources with 
/ > 10~ 8 ph/(s cm 2 ), i.e. just an order of magni- 
tude weaker than the EGRET limit. Thus, with the 
sensitivity improved by an order of magnitude be- 
low the EGRET threshold, both models would be 



5. CONCLUSIONS 



Simple, discrete source models of the GRB have been 
constructed to investigate effects of the background 
anisotropies at different angular scales. Source 
counts predicted in the balazar and starburst galaxy 
models are distinctly different. We have shown that 
precise measurements of the GRB using of the IN- 
TEGRAL detectors would allow to distinguish be- 
tween these two cases. Isotropic GRB would indicate 
that the large number of faint sources contributes to 
the GRB making the starburst model more likely, 
while detectable GRB fluctuations would point to 
the blazar model. 
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